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A new  technique  to  reduce  the  effect  of  quantization  in  PCM  image 
coding  is  presented  in  this  report.  The  new  technique  consists  of  Roberts' 
pseudonoise  technique  followed  by  a noise  reduction  system.  The  technique 
by  Roberts  effectively  transforms  the  signal  dependent  quantization  noise 
to  a signal  Independent  additive  random  noise.  The  noise  reduction  system 
that  follows  reduces  the  additive  random  noise.  Some  examples  are  given 
to  Illustrate  the  performance  of  the  new  quantization  noise  reduction 
system. 
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I.  INTRODUCTION 


In  coding  an  Image  by  a Pulse  Code  Modulation  (PCM)  technique,  an  image 
is  sampled  and  a bit  rate  reduction  is  achieved  by  reducing  the  number  of 
quantization  levels  assigned  to  each  sample.  If  the  number  of  quantization 
levels  is  not  sufficiently  large,  noticeable  gray  scale  contouring  results 
due  to  a luminance  jump  in  the  reconstructed  image  in  a region  where  the 
original  image  luminance  varies  slowly.  For  typical  monochrome  Images, 
experience  (1,2)  has  showti  that  at  least  five  bits  per  sample  are  generally 
necessary  to  avoid  gray  scale  contouring.  If  the  degradation  due  to  quan- 
tization can  be  reduced,  then  Images  of  equivalent  quality  may  be  generated 
with  a smaller  number  of  quantization  levels.  Thus,  it  is  of  Interest  to 
develop  techniques  to  reduce  the  degradation  due  to  quantization. 

Various  techniques  have  been  proposed  (1,3, A, 5)  in  the  literature  to 
reduce  the  visual  effect  of  gray  scale  contouring.  In  most  cases,  the 
general  approach  has  been  to  transform  gray  scaling  contouring  to  a 
different  type  of  degradation  which  Is  more  easily  tolerated  by  human 
viewers.  One  such  example  is  Roberts'  pseudonoise  technique  (3).  Even 
though  such  an  approach  has  been  successful  In  eliminating  gray  scale  con- 
touring, the  quantization  noise  In  the  mean  square  sense  generally  Increases 
and  the  degradation  due  to  quantization  may  be  less  annoying  but  neverthe- 
less Is  quite  visible.  In  this  paper,  we  develop  a new  quantization  noise 
reduction  technique  by  attempting  to  eliminate  the  quantization  noise. 

Gray  scale  contouring  due  to  quantization  is  a signal  dependent  de- 
gradation and  therefore  it  is  difficult  to  apply  signal  processing  techniques 
for  its  elimination.  Roberts  (3),  however,  has  shown  that  gray  scale  con- 
touring due  to  quantization  can  be  transformed  to  a signal  independent 
additive  random  noise  with  a slightly  larger  mean  square  error  by  adding  a 
known  random  noise  to  the  sampled  image  before  quantization  and  subtracting 
the  same  random  noise  after  quantization.  Once  the  quantization  effect  is 
converted  to  a signal  independent  additive  random  noise,  then  the  quanti- 
zation noise  may  be  more  easily  dealt  with  by  currently  available  signal 


1 


processing  techniques.  This  Is  the  basis  for  the  new  quantization  noise 
reduction  technique  presented  In  this  paper.  As  will  be  seen  later,  such 
an  approach  to  reducing  the  quantization  noise  has  the  potential  to  elim- 
inate gray  scale  contouring,  reduce  the  mean  square  error,  and  generate 
higher  quality  Images  than  existing  techniques. 

The  overall  objectives  of  this  paper  are  to  present  a new  technique 
to  reduce  the  quantization  noise  In  PCM  Image  coding  and  Illustrate  Its 
performance.  In  Section  II,  we  describe  a basic  PCM  Image  coding  technique, 
Roberts'  pseudo-noise  technique  and  a new  quantization  noise  reduction 
technique  that  have  been  Implemented  for  their  performance  comparison.  In 
Section  III,  we  compare  and  discuss  the  performance  of  the  techniques 
described  In  Section  II. 


II.  A NEW  TECHNIQUE  FOR  QUANTIZATION  NOISE  REDUCTION 


In  this  section,  we  develop  a new  technique  to  reduce  the  quantization 
noise  in  Pt^  image  coding.  We  also  describe  a basic  Ptn  coding  technique 
and  Roberts'  pseudonoise  technique  which  will  be  compared  in  performance 


with  the  technique  developed  In  this  section. 
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The  basic  POI  image  coding  technique  that  has  been  implemented  is  showi\ 
in  Figure  1.  In  the  Figure,  f(nj,n2)  represents  a noise-free  digital  image 
and  q(nj,n2)  represents  an  image  quantized  by  a uniform  quantizer. 


UNIFORM 
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Fig.  1.  A Basic  PilM  Coding  Technique 


In  the  pseudonoise  technique  by  Roberts,  a known  random  noise  is  added 
to  the  noise-free  image  f(nj,n,)  before  quantization  and  the  same  ramh'm 
noise  Is  subtracted  after  quantization,  as  is  shown  in  Figure  2.  In  the 
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Roberts'  Pseudonolse  Technique 


ftt^ure,  w(n^,n^)  represents  white  noise  generated  by  a uniform  probability 


density; 
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where  A represents  the  spacing  between  reconstruction  levels.  The  sequence 
r(nj^tn^)  represents  the  reconstructed  image. 

The  new  technique  for  quantization  noise  reduction  consists  of  Roberts' 
pseudonoise  technique  followed  by  a noise  reduction  system,  as  is  showti  in 
Figure  3.  The  sequence  s(nj^,n2)  represents  the  reconstructed  image.  The 
technique  by  Roberts  effectively  transforms  the  signal  dependent  quantiza- 
tion noise  to  a signal  independent  additive  random  noise  (3).  The  noise 
reduction  system  that  follows  reduces  the  additive  random  noise.  For 
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Fig.  3.  A New  Technique  for  Quantization  Noise  Reduction 
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the  noise  reduction  system,  we  have  Implemented  on  a short  space  basis  the 
"Spectral  Subtraction  Image  Restoration"  (SSIR)  technique  recently  proposed 
by  Llm  (6).  The  SSIR  technique  was  chosen  due  to  its  effectiveness  in 
reducing  additive  random  noise  without  significant  signal  distortion. 

In  the  short  space  SSIR  technique,  the  image  r(nj,n2)  in  Figure  3 is 
divided  into  many  subimages  each  of  which  is  restored  separately  and  then 
the  restored  Images  are  combined  to  form  s(nj^,n2).  More  specifically,  let 
rCOj^.n^)  be  represented  by 


r(nj^,n2)  = fCn^^.n^)  + d(nj^,n2) 


(2) 


where  d(nj^,n2)  denotes  the  additive  random  noise  that  results  from  Roberts' 
pseudonoise  technique.  By  applying  a 2-D  window  function  w^j(n^,n2)  to 
equation  (2), 

r(n^,n2)  • Wjj(nj,n2)  » f(nj^,n2)  • w^j(nj^,n2)  + d(nj^,n2)  • w^^(n^,n2)  • (3) 
Rewriting  equation  (3), 


(4) 


where  represents  r(nj^,n2)  • Wj^j(n^,n2),  and 

*^ij^*'l’"2^  are  similarly  defined.  To  estimate  the  noise-free  subimage 
^lj^"l’"2^  from  in  equation  (4),  F^j  (u)j^,a)2) , the  discrete  space 

Fourier  transform*  of  f^j(nj,n2),  is  first  estimated  and  then  inverse  Fourier 

* The  definition  of  discrete  space  Fourier  transform,  power  spectrum 
and  energy  spectrum,  and  the  determination  of  the  normalization 
constant  "k"  can  be  found  in  references  (6)  and  (7). 
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transformed.  The  discrete  space  Fourier  transform  is  estimated  by 

a particular  form  of  spectral  subtraction  (6); 


’'lj^‘*’r‘^2^  ■ (|Ry(Wj,W2)|  - . e-^  <Rjj(Wj,u-2) 


£ 

for 


and  0 othervlse 


where  represents  an  estimate  of  F^j  (u)^ ,u),) , represents 

the  discrete  space  Fourier  Transform  of  r^^Cn^.n^),  < represents 

the  phase  of  R^j  (a)j^,u)2) , "oi"  is  a constant,  "k"  Is  a scaling  factor  that 

normalizes  the  power  and  energy  spectral  densities,  and  P.iw, .w-)  represents 

u 1 2 

the  power  spectrum  of  the  additive  random  noise.  From  the  estimated 
^ij  ’*^2^  ’ Figure  3 Is  obtained  by  combining  the  restored 

subimages; 


sCn^^.n^) 


EE 

1 j J 


where  represents  the  estimated  f^j(nj^,n2). 

In  implementing  the  short  space  SSIR  technique  dlscvissed  above,  a 
separable  2-D  triangular  window  of  size  16  x 16  pixels  overlapped  with  its 
neighboring  window  by  half  the  window  duration  In  each  dimension  was  used 
for  w^j(nj^,n2)  value  of  "a"  was  assumed  to  be  2.  The  additive 

random  noise  d(nj^,n2)  was  approximated  as  white  noise  with  the  same 


j 


statlstios  as  w(n  ,n^)  In  equation  (1),  and  thus  the  power  spectrum  T (a)  ,u) 

*■  2 d I 

was  approximated  to  be  flat  with  t tie  spectral  amplitude  of  A /12  wliere  A 

represents  tlie  spacing  between  reconstruction  levels.  Furtlier  details  on 

the  theoretical  development.  Implementation  and  performance  of  ttie  SSIK 

technique  can  be  found  in  (b). 
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In  this  section,  we  present  a few  examples  to  illustrate  the  performance 

I 

[•  of  the  quantization  noise  reduction  technique  developed  In  Section  II.  In 

i Figure  4 are  shown  three  original  Images  of  256x2‘j6  pixels  with  each  pixel 

represented  bv  8 bits.  The  three  Images  will  be  referred  to  as  Image  1,2 

i 

and  3.  In  Figure  5 are  shown  Che  reconstructed  images  from  Image  1. 
j Specifically,  in  Figures  5(a),  (b)  and  (c)  are  sliown  the  reconstructed  images 

> q(nj,n.,)  by  the  basic  PCM  coding  technique  in  Figure  1,  r(nj,n.,)  by  Roberts' 

[ pseudonoise  technique  in  Figure  2 and  3(11^11.,)  by  tlie  quantization  noise 

[ reduction  technique  in  Figure  3,  with  three  bits  per  s.imple  assigned  to  the 

uniform  quantizer.  Figures  5(d),  (e)  and  (f)  are  equivalent  to  Figures  5(a), 
(b)  and  (c)  except  that  two  bits  per  sample  were  assigned  to  the  uniform 
quantizer.  Figures  6 and  7 are  equivalent  to  Figure  5 except  that  the 
Images  were  reconstructed  from  Images  2 avtd  3.  respectively.  From  the  above 
examples.  It  is  clear  that  Roberts'  pseudoni'ise  technique  eliminates  gray 
scale  contouring  but  the  reconstructed  images  appear  to  be  noisy.  Ky 
cascading  an  effective  noise  reduction  system,  however,  noticeable  noise 
reduction  can  be  achieved  without  significant  signal  distortion. 

Even  though  a lower  moan  square  error  does  not  necessarily  imply  lilgher 
image  quality,  the  mean  square  error  is  a useful  measure  in  that  it  is  easy 
to  compute  and  has  at  least  some  correlation  with  image  quality.  For  each 
of  the  examples  Illustrated  in  this  section,  a normalized  mean  square  error 
(NMSE)  was  computed  and  is  tabulated  in  Table  1.  The  MMSE  in  between  an 
original  image  f(nj^,n.,)  and  a reconstructed  image  fin^.n^)  is  defined  by 

))  (f  (n^,n.,)-f  ' (n  j ,n.,)  ) “ 

"1  "2  , 

NM.SE(f(n  ,n.,);f(n  ,n^))  = 100  

E T (f  (iij  .n.^)-El  t (iij  ,!!.,)])*■ 


Fig.  5.  (a)  3 bit  PCM  System 


Cd)  2 bit  PCM  System 


(b)  3 bit  Roberts'  Pseudonoise  System  (e)  2 bit  Roberts'  Pseudonoise  System 


(c)  3 bit  New  Quantization  Noise 
Reduction  System 


(f)  2 bit  New  Quantization  Noise 
Reduction  System 
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(c)  3 bit  New  Quantization  Noise 
Reduction  System 


(f^  2 bit  New  Quantization  Noise 

Reduction  System 


Fig.  6.  (a)  3 bit  PCM  System 


(d)  2 bit  PCM  System 


(b)  3 bit  Roberts'  Pseudonoise  System  (e)  2 bit  Roberts'  Pseudonoise  System 


Fig.  7.  (a)  3 bit  PCM  System 


(b)  3 bit  Roberts'  Pseudonoise  System  (e)  2 bit  Roberts'  Pseudonoise  System 


(f)  2 bit  New  Quantization  Noise 
Reduction  System 


(c)  3 bit  New  Quantization  Noise 
Reduction  System 


whort*  f’Crij.n^)  ■ a.f(n|,n^)+b  with  "a"  and  "h"  I'hosen  sui'h  that 
Klf(n  ,n^)l  - K|f'(n  ,n,)I  and  Var  (fin  ,n  ))  • Var  If'(n  ,n  )).  K|f(n  ,n  )] 
and  Var  (fln^.n^)!  rt'presont  nK>an  and  variance  of  f(nj,n^),  and  Kl  f ’ (n^  ,n^> ) 
and  Var  If'ln^.n^))  are  Himllarly  defined.  The  definition  of  NMSK  Klven  hv 
equation  ( ))  has  the  property  that  Ilnearlv  sealing  and  adding  a hlas  to 
f(n  .n  ) or  f(n  ,n  ) do  not  affect  NM.SK  (f(n,,n.);  f(n,,nj),  and  NMSK 
(f(nj,n,);  ffn^.n^))  equals  NMSK  ( f (n  j . n^) . Kn^.n^).  Kri>m  Table  1.  (t  Is  clear 
that  Images  reconstructed  hv  the  quantization  noise  reduction  technique 
developed  In  this  pape>"  have  significantly  lower  NMSK  than  Images  re- 
ci'nstructed  by  a basic  P(>1  coding  technique  or  Roberts'  pseudonoise  teclinl<tue. 

In  this  paper,  we  have  considered  reducing  the  quantization  noise  bv 
first  transforming  gray  scale  contouring  to  .an  additive  random  noise  .and 
then  applying  a noise  reduction  system  to  reduce  the  random  noise.  Kv«'n 
though  a detailed  subjective  evaluation  has  nvit  been  performed  to  measure 
the  inuige  qualltv  Improvement,  the  results  In  this  paper  are  quite  en- 
couraging and  Indicate  that  such  an  approach  to  reduce  the  quantization 
noise  In  PCM  Image  coding  can  potentially  aid  In  developing  new  image  coding 
techniques  or  Improving  existing  coding  techniques  (8,9)  which  Incorpvirate 
PCM  coding  as  an  Integral  part. 
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